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1. Introduction 



The expression for the magnetic flux density from a 
single rectangular loop of wire of many turns can be 
found in text books and various publications [1-3]. The 
rectangular geometry is convenient, in part, because the 
expressions for the three spatial components of the flux 
density are in closed form. Single square coils have been 
used for calibration of extremely low frequency mag- 
netic field meters for applications that require uncertain- 
ties of a few percent [2]. Multiple rectangular loops with 
a common axis have found applications in a number of 
fields, including biological exposure systems for in vivo 
and in vitro studies [3,4]. It is also noteworthy that a 
square Helmholtz coil produces a greater volume of 
nearly uniform magnetic field than a circular Helmholtz 
coil of comparable dimensions [5]. This paper develops 
expressions for the magnetic flux density produced by 
three rectangular loops of wire that lie in the same plane, 
i.e., loops that are not co-axial. The geometry is similar 
to that used in some walk-through metal detectors and 



electronic article surveillance systems. By inspection, 
the expressions for more or fewer loops are easily deter- 
mined. We consider static and time varying fields that 
are quasi- static. In the latter case, the wavelength A of 
the time varying field is much greater than any dimen- 
sion or distance of interest. For example, a 1 MHz alter- 
nating field (A ~ 300 m) is well approximated as being 
quasi- static a few meters or less from loops of compara- 
ble dimensions. The quasi-static condition allows us to 
solve the static field problem first and, with negligible 
error, introduce the time dependence as a multiplicative 
factor, e.g., the direct current in the field equations could 
be replaced with an alternating current. The field equa- 
tions are for rectangular loops with a single turn of wire. 
The magnetic flux density for loops with more than one 
turn are found by multiplying the equations by the ap- 
propriate number of turns. 
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2. Field Equations 

We follow the development of Weber [1] by first con- 
sidering the vector potential for a rectangular loop of 
wire in the x-y plane, A^ and Ay , and then calculating the 
vector components of the magnetic flux density using 
the relations 



The z -component of the magnetic flux density at 
P{x,y,z) is 



f^ply 
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For a single rectangular loop of wire of negligible wire 
cross section, designated as loop 1, with side dimensions 
2(21 by 2bi as shown in Fig. 1, the components of the 
vector potential are [1] 
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where /j,o is the magnetic constant (also called the mag- 
netic permeability of vacuum), and I\ is the current in 
the loop. 

The parameters r^ r2, r^, and r^, are the distances from 
the corners of the loop to the point P(x ,y,z) where the 
magnetic flux density will be evaluated (see below and 
Fig. 1). 



P(x,y,z) 




Fig. 1. Geometry for a single rectangular loop of wire in the x-y 
plane. The magnetic flux density is evaluated at point V(x, y,z)- 



Equation (4) is equivalent to that given in Ref. [1], but 
perhaps in a more convenient form for writing a com- 
puter program to calculate the magnetic flux density. 

From Eqs. (1) to (3), the expressions for the x- and 
y -components of the magnetic flux density can be read- 
ily derived and are 



and 
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The X -component of the vector potential for a second 
loop of wire of side dimensions 2a2 by 2/^2 that is dis- 
placed from the origin by a distance 52 and bisected by 
the y-axis (see Fig. 2) is given by [1] 
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and I2 is the current in loop 2. 

The integrals can be solved using elementary methods 
and yield 
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,P(x,y,z) 




Fig. 2. Geometry for a second rectangular loop of wire in x-y plane. 
The point P(x, y,z) coincides with that in Fig. 1 (note that the scales 
of Figs. 1 and 2 are not the same). 






(r{ + 02+ x) (ri — 02 + x) 



(r{ — a2 + x) (rl + a2 + x) 



(8) 



where r{ , ri , ri , and ri are the distances from the corners 
of loop 2 to the point P{x,y,z) where the magnetic flux 
density will be evaluated (see below). 

The expression for Ay2 can be similarly determined 
and is given by 
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Taking the appropriate derivatives of Eqs. (8) and (9), 
the expression for the z -component of the magnetic flux 
density at F(x,y,z) associated with loop 2 is 
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From Eqs. (1), (8), and (9), the x- and y -components of 
the magnetic flux density due to loop 2 are 
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The equations for the flux density components at 
P(x ,y,z) from a third rectangular loop with side dimen- 
sions 2fl3 by 2i>3, displaced from the origin by a distance 
S3 and bisected by the y -axis follow by inspection. That 
is 
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and h is the current in loop 3. 

The X- and y -components of the magnetic flux den- 
sity due to loop 3 are 
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The spatial components of the magnetic flux density 
at ¥{x,y,z) due to all three loops (Fig. 3) are found by 
summing the respective contributions from each loop, 
i.e.. 
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P(x,y,z) 



Fig. 3. Geometry for three rectangular loops of wire in a vertical 
plane. The origin of the coordinate system is at the center of loop 1 . 
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For direct currents in the loops, the direction of the 
magnetic flux density will remain fixed and is de- 
scribed by the vector 



B=B,Ti + ByTJ + B,Tk, 



(17) 



where i,j, and k are unit vectors along the x, y, and z 
directions, respectively. The magnitude of the magnetic 
flux density vector will also be constant and equal to 



\B\ = \/B}r + B^T + Bi 



(18) 



For alternating currents in the loops that are in phase, 
for example /isin(wf), /2sin((yf)j and I^&m{(i)t), the 



magnetic flux density is described by the vector 

B = [B,T i + B,TJ + B,Tk]smicot). (19) 

where /i, h, and I3 are current amplitudes, co is the 
angular frequency, and t is the time. The flux density is 
said to be linearly polarized because of its oscillatory 
motion along a straight line. The magnitude of the vec- 
tor will be time dependent and equal to 



iVfi 



^T + Byl + Bjr sin(wf) 



If the alternating currents in the various loops are not 
in phase, the magnetic flux density vector will rotate 
and the point of the vector will, in general, trace an 
ellipse [6] . The magnitude and direction of the magnetic 
flux density at a given point in space will change as a 
function of time. For this case, the flux density is said 
to be elliptically polarized. 

As a convenience to the reader, a program for calcu- 
lating the static magnetic flux density from three coils 
in the x-y plane as shown in Fig. 3 is provided in Ap- 
pendix A. 



3. Appendix A. Program for Calibrating 
the Magnetic Flux Density 

The Quick Basic' program below calculates the mag- 
netic flux density at a point x,y,z specified in the input 
statement for three rectangular loops of single turn wire 
with side dimensions 2ai, 2bi, 2a2, 2i>2, 2fl3, and Ib^. The 
loops are separated by distances Sz and S3, and are bi- 
sected by the y -axis. The origin of the coordinate system 
is at the center of loop 1 . Figure 3 shows the geometry 
of the three loops and a point P(x, y, z) where the 
magnetic flux density is calculated. The program as- 
sumes that Ii= l2 = h = I Aof direct current and that the 
current is in the counter clockwise direction as one views 
the current looking down on the x-y plane (e.g., see Fig. 
1). All dimensions are in meters and the magnetic flux 
density will be in units of 10"' T. For example, with la,, 
2bi, 2a2, 2b2, 2^3, and 2^3 equal to 1.0 m, 0.5 m, 2.0 m, 
1 m, 0.5 m, and 0.3 m, respectively, separated by dis- 
tances .?2 = 1-5 m and S3 = 2.5 m, the magnitude of the 
magnetic flux density IB I at x = 0.75 m, y = 0.75 m, and 
z= 1 mis 0.95170 X 10"' T. 



Certain commercial products are identified in this paper to foster 
understanding. Such identification does not imply recommendation or 
endorsement by the National Institute of Standards and Technology, 
nor does it imply that the materials identified are necessarily the best 
available for the purpose. 
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REM PROGRAM TO CALCULATE B-FIELD FROM THREE RECTANGULAR SINGE TURN COILS WITH 
REM SIDE DIMENSIONS 2AI X 2B1, 2A2 X 2B2, AND 2A3 X 2B3, AND SEPARATED BY DISTANCES 
REM S2 AND S3, RESPECTIVELY. ALL THREE LOOPS ARE IN THE X-Y PLANE AND THE 
REM Y-AXIS BISECTS EACH RECTANGLE; SEE FIG. 3. 

DEFDBL A-H, K-Z 
CLS 

DIM C(4), D(4), P(4), Q(4), T(4), R(4) 

DIM CP(4), DP(4), PP(4),QP(4), TP(4), RP(4) 

DIM CPP(4), DPP(4), PPP(4), QPP(4), TPP(4), RPP(4) 

PRINT "ENTER: AI, Bl, A2, B2, A3, B3, S2, S3, X, Y, Z" 
INPUT AI, Bl, A2, B2, A3, B3, S2, S3, X, Y, Z 

REM TERMS FOR LOOP 1 

C(l) = AI +X 

C(2) = AI - X 

C(3) =-C(2) 

C(4) =-C(l) 

D(I) = B1 + Y 

D(2) = D(I) 

D(3) = Y - Bl 

D(4) = D(3) 

R(l) = SQR(C(I)'^2 + D(I)'^2 + Z'^2) 

R(2) = SQR(C(2)'^2 + D(I)'^2 + Z'^2) 

R(3) = SQR(C(2)'^2 + (Bl- Y)'^2 + Z'^2) 

R(4) = SQR(C(1)'^2 + (Bl- Y)'^2 + Z'^2) 

REM CORRESPONDING TERMS FOR LOOP 2 

CP(1) = A2 + X 

CP(2) = A2 - X 

CP(3) =-CP(2) 

CP(4) =-CP(l) 

DP(1) = Y + B2 - S2 

DP(2) = DP(1) 

DP(3) = Y - B2 - S2 

DP(4) = DP(3) 

RP(1) = SQR(CP(1)'^2 + DP(ir2 + Z'^2) 

RP(2) = SQR(CP(2)'^2 + DP(1)'^2 + Z'^2) 

RP(3) = SQR(CP(2)'^2 + DP(3)'^2 + Z'^2) 

RP(4) = SQR(CP(1)'^2 + DP(4)'^2 + Z'^2) 

REM CORRESPONDING TERMS FOR LOOP 3 

CPP(l) = A3 + X 

CPP(2) = A3 - X 

CPP(3) =-CPP(2) 

CPP(4) =-CPP(l) 

DPP(l) = Y + B3 - S3 

DPP(2) = DPP(l) 

DPP(3) = Y - B3 - S3 

DPP(4) = DPP(3) 

RPP(l) = SQR(CPP(1)'^2 + DPP(1)'^2 + Z'^2) 

RPP(2) = SQR(CPP(2)'^2 + DPP(1)'^2 + Z'^2) 

RPP(3) = SQR(CPP(2)'^2 + DPP(3)'^2 + Z'^2) 

RPP(4) = SQR(CPP(1)'^2 + DPP(4)'^2 + Z'^2) 
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REM CALCULATE X COMPONENT OF LOOP 1 (BXl) 

FOR J = 1 TO 4 

P(J) = (Z*(-1)'^(J+1))/((R(J) + D(J))*R(J)) 

NEXT J 

BXl = 

FOR I = 1 TO 4 

BXl = BXl + P(I) 

NEXT I 

PRINT "BXFIELDl ="; BXl 

REM CALCULATE Y COMPONENT OF LOOP 1 (BYl) 

FOR J = 1 TO 4 

P(J) = (Z*(-ir(J+l))/((R(J) + ((-l)'^(J +1))*C(J))*R(J)) 

NEXT J 

BYl = 

FOR I = 1 TO 4 

BYl = BYl + P(I) 

NEXT I 

PRINT "BYFIELDl ="; BYl 

REM CALCULATE Z COMPONENT OF LOOP 1 (BZl) 

FOR J = 1 TO 4 

P(J) = ((-irJ)*D(J)/(R(J)*(R(J) + ((-l)'^(J + 1))*C(J))) 

Q(J) =-C(J)/(R(J)*(R(J) + D(J))) 

T(J) = P(J) + Q(J) 

NEXT J 

BZl =0 

FOR I = 1 TO 4 

BZl = BZl + T(I) 

NEXT I 

PRINT "BZFIELDl ="; BZl 

PRINT 

REM CORRESPONDING TERMS FOR LOOP 2; BEGIN WITH X COMPONENT OF LOOP 2 (BX2) 

FOR J = 1 TO 4 

PP(J) = (Z*(-1)'^(J+1))/((RP(J) + DP(J))*RP(J)) 

NEXT J 

BX2 = 

FOR I = 1 TO 4 

BX2 = BX2 + PP(I) 

NEXT I 

PRINT "BXFIELD2 ="; BX2 

REM CALCULATE Y COMPONENT OF LOOP 2 (BY2) 

FOR J = 1 TO 4 

PP(J) = (Z*(-1)'^(J+1))/((RP(J) + ((-l)'^(J +1))*CP(J))*RP(J)) 

NEXT J 

BY2 = 

FOR I = 1 TO 4 

BY2 = BY2 + PP(I) 

NEXT I 

PRINT "BYFIELD2 ="; BY2 
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REM CALCULATE Z COMPONENT OF LOOP 2 (BZ2) 

FOR J = 1 TO 4 

PP(J) = ((-1)AJ)*DP(J)/(RP(J)*(RP(J) + ((-l)'^(J + 1))*CP(J))) 

QP(J) =-CP(J)/(RP(J)*(RP(J) + DP(J))) 

TP(J) = PP(J) + QP(J) 

NEXT J 

BZ2 = 

FOR I = 1 TO 4 

BZ2 = BZ2 + TP(I) 

NEXT I 

PRINT "BZFIELD2 ="; BZ2 

PRINT 

REM CORRESPONDING TERMS FOR LOOP 3; BEGIN WITH X COMPONENT OF LOOP 3 (BX3) 

FOR J = 1 TO 4 

PPP(J) = (Z*(-1)'^(J+1))/((RPP(J) + DPP(J))*RPP(J)) 

NEXT J 

BX3 =0 

FOR I = 1 TO 4 

BX3 = BX3 + PPP(I) 

NEXT I 

PRINT "BXFIELD3 ="; BX3 

REM CALCULATE Y COMPONENT OF LOOP 3 (BY3) 

FOR J = 1 TO 4 

PPP(J) = (Z*(-1)'^(J+1))/((RPP(J) + ((-ir(J +1))*CPP(J))*RPP(J)) 

NEXT J 

BY3 = 

FOR I = 1 TO 4 

BY3 = BY3 + PPP(I) 

NEXT I 

PRINT "BYFIELD3 ="; BY3 

REM CALCULATE Z COMPONENT OF LOOP 3 (BZ3) 

FOR J = 1 TO 4 

PPP(J) = ((-1)'^J)*DPP(J)/(RPP(J)*(RPP(J) + ((-l)'^(J + 1))*CPP(J))) 

QPP(J) =-CPP(J)/(RPP(J)*(RPP(J) + DPP(J))) 

TPP(J) = PPP(J) + QPP(J) 

NEXT J 

BZ3 = 

FOR I = 1 TO 4 

BZ3 = BZ3 + TPP(I) 

NEXT I 

PRINT "BZFIELD3 ="; BZ3 

PRINT 

REM COMBINE SPATIAL COMPONENTS OF B-FIELD FROM LOOPS I, 2, AND 3 

BXT = BXl + BX2 + BX3 

BYT = BYl + BY2 + BY3 

BZT = BZl + BZ2 + BZ3 

REM MAGNITUDE OF B VECTOR 

BMAG = SQR(BXT'^2 + BYT'^2 + BZT'^2) 

PRINT "BXTOTAL ="; BXT; "BYTOTAL ="; BYT; "BZTOTAL ="; BZT 

PRINT "MAGNITUDE OF B ="; BMAG 

PRINT "X ="; X; "Y ="; Y; "Z ="; Z 

END 
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